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A variational approach is used to study the superconductor-insulator transition in two-band gran- 
ular superconducting films using a resistance-shunted Josephson junction array model in this letter. 
We show that a zero-temperature metallic phase may exist between the superconducting and insu- 
lator phases which is absent in normal single band granular superconducting films. The metallic 
phase may be observable in some dirty pnictide superconductor films. 
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Intensive studies had been devoted to the problem of 
superconductor-insulator (SI) transition in low- Tc thin 
films. These systems undergo phase transitions from su- 
perconductor to insulator as a function of disorder, film 
thickness as well as external magnetic fields |T]. The SI 
transition is usually modelled by a Josephson junction 
array model, expressed in terms of the phases of the su- 
perconductor order parameter 9i's on superconducting 
grain i's. The Hamiltonian describing the system con- 
sists of the Josephson coupling between superconducting 
grains J cos{9i — 9j) where (i.j) are nearest neighbor 
sites, and the charging energy ~ ^{di — ^j)^- The sys- 
tem is in a superconducting phase if the Josephson term 
dominates, and is in the insulator phase if the charg- 
ing energy dominates. It has been proposed by different 
authors that a dissipative term arising from coupling be- 
tween superconducting grains and a dissipative metallic 
bath may also be important in describing the SI transi- 
tion (shunted Josephson array model) ^-tdj . In particular, 
a zero-temperature metallic phase between superconduc- 
tor and insulator phases may be stabilized by dissipation. 

The physical reason behind the metallic phase is as 
follows: Imagine first a state dominated by charging en- 
ergy. In this case the metallic bath would screen the 
Coulomb potential, leading to a weakening of charging 
energy and drives the system towards a metallic phase 
if the resistance is small enough {R < i?c/)0- Alter- 
natively, the coupling of Cooper pairs in the supercon- 
ducting phase to a dissipative environment suppresses 
coherent tunnelling of Cooper pairs between grains ow- 
ing to the Calderia-Leggett effect [B] and superconduct- 
ing coherence is destroyed if i? > RcS- As a result a 
metallic phase between the superconducting and insulat- 
ing phases may exist if Rcs < R < Rd- The metal- 
lic phase, if exist, is a new phase of matter because of 
participation of incoherent boson (Cooper pairs) in low 
temperature transports which is absent in usual metals. 
Experimentally the zero-temperature metallic phase in 
single band superconducting films has not been found to 
exist so far in the absence of external magnetic fields, 
consistent with a theoretical finding that Rd < RcS hi 
single-band superconductors [7]. 



More recently, superconductors with more than one or- 
der parameters, i.e., the multi-band superconductors [8] 
have raised attention in the physics community. Exam- 
ples of multi-band superconductors include MgB2[9j and 
the pnictide superconductors [10]. It is interesting to see 
whether a metallic phase may exist more easily between 
the Si-transition in these materials. This is the purpose 
of this letter. 

Using a variational approach, we consider in this let- 
ter the superconductor-insulator transition in two-band 
s (and s-|-)-wave superconducting films where the possi- 
bility of an intermediate metallic phase is investigated. 
We show that contrary to the case of single-band su- 
perconductors, a physically realizable condition for the 
zero-temperature intermediate metallic phase is found 
for these systems. We propose that the metallic phase 
may be observable in some recently discovered disordered 
pnictide superconductors[ll . 
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Figure 1: A qualitative sketch of our two band Josephson 
junction model. The two bands on different grains are con- 
nected by a capacitor, a resistor and Josephson coupling (not 
shown in the sketch). The two bands on the same grain are 
connected by in-grain inter-band Josephson coupling Ji 

We start with the phase-action which is a generaliza- 
tion of the phase action used to study superconductor- 
insulator transition in one-band systems [31 |31 [7]. The 
system is schematically sketched in Fig.l. The action 
describes a resistance-shunted Josephson network of two- 
band superconductor grains and is given in imaginary 
time hy S = Sg + S'diss: where 
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5e = V r drilCabiAjff - JabCosiA^ef)] 

, "'0 



is the phase action without the dissipative term. a,b — 
1, 2 represent the two different bands in a grain, and Of 
is the phase of band a superconducting order parame- 



ter in grain i. A^Of^ = 0f — 0^^^ represents the phase 
difference between band a and b superconducting order 
parameters in neighboring grains i, i + v, respectively 
and Jab > is the corresponding Josephson couphng en- 
ergy. ^CabiAi.Of'')'^ represents the charging energy aris- 
ing from charge imbalance between band a and band b 
electrons on grain i and i+v respectively, where Cab is the 
corresponding capacitance. Jj is the in-grain inter-band 
Josephson coupling which favors 0f = 9l + tt for J/ > 0, 
leading to a s± superconductor and favors 9} = 9f for 
J/ < (s-wave superconductor). 



Sdi: 



J: dr dr' 



xaab(T- T )sin 



2Q 



(2) 



where a°b(r) = (/i/4e2i?°J[T/sin(7rTr)]2. is the re- 
sistance between band a and band b electrons on grains 
i and i + i^, respectively (see Fig.(l)) and Q = 2 is the 
charge of a Cooper pair. 0,^0}^ = D^Of^ = [A^6'," -I- 
A^9f% D^Of = [|A^6'r-f iA^6'f ], where 1(2) = 2(1). 
Sdiss is derived phenomenologically from a multi-band 
resistance network model represented by Fig.l. The de- 
tails of the derivation can be found in the supplementary 
materials. 

To simplify calculation we shall consider the grains 
forming a two-dimensional square lattice with J12 = 
in our following analysis. With the later condition the 
s and s± superconductors can be transformed to each 
other by simply shifting 0^ 9f + it. The main effect 
of J12 is to renormalize Jj ^ Jj — zJi2 where z is the 
lattice co-ordination number and is not going to affect 
our conclusion in renormalization-group sense. 

Due to the compactness of the phase field (e*(^+2n7r) _ 
e*^), the phase variables (r) can be decomposed into a 
periodic part and a winding number contribution, 



27rnfr 



where 0fo(/^) = ^io(O) ^^"^ ''^f '^^^ ^'^y arbitrary integer 
(winding number). With this decomposition the phase 



action becomes 
27r2 



Se ^ 



Y: CabA.nf +T.^J d-i^^^tof 

i.u,a,h i.u.a.h ^ 

^Y.J-- dT cos[A,9fS + -yA^nT] 

i,i/,a ^ 

2'KT 

+JiY^ dTy^cos[A9,o + —An,]. (3) 
. "'0 „ P 



where A^nf = nf - A9,o = 0}o{t) - 91{t) and 

An, — n} — and 



'S'diss ~^ 



a°,(r-r')cos[?^^i^^i?.nf] 
x[D,9t,\T)-D^9'^^{T')r. 



(4) 



where D^, nf is defined in the same way as D^9f with 
Aj,6'f^ — >■ A^nf = nf — n^^^^. We have assumed strong 
dissipation and keep only to second order terms of Ai^flgf 
in S'diss for simplicity [Tj- 

To proceed further we employ a varia- 
tional approach T^. We consider a trial action 
•Striai — '5't^iai + '^'triai whcre the pcriodic and the 
winding number contributions to 9 are decoupled. 



^t^iai = E fdrffiAJt) 

2, 1/, a, 6 ^ 

7-cff r!3 

+ dr{A9,o) 

- y 



Tcff 

2 I ''aa 



/3 /-/J 

dr / dT'af^°(T-r') x 

nab ( ^\ T~) nab /'^f\'\2 



[D^9t^ir)-D,9tS{r')Y 



(5) 



is an effective action describing Gaussian fluctuations of 
the periodic phases around the saddle point 91q{t) — 0, 
6»2g(r) = 0(7r) for J/ < (>)0 and 
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trial 



2^ 



E a.Ax" -/3Jf^E^K'-"') 

i,u^a^b i 

Qtt 



4i?' 



ab I 



(6) 



in an effective action for the winding number field. S'^J.jj^j 
is a generalized absolute solid-on-solid model (ASOS) for 
two species of winding numbers with additional /3J^'^^ 
terms originating from superconductivity. J^^^Jf^, Jaa^ i 
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Jj^^, af^°{T)/a{T) = R°a/Rf'° are variational param- 
eters to be determined by minimizing the free energy of 
the system given approximately by F = Fq + (S" — S'triai), 
where Fq is the free energy computed using S'tiiai and (...) 
denotes averages taken with respect to S'triai- 

The different phases can be identified in our trial ac- 
tion as follows: First we note that 5'^;^! describes a stable 
superconducting phase as long as the phase stifFness's 
satisfy J^^ > 0. The nature of the J^^ = (non- 
superconducting) state is determined by •S'^^jai which de- 
scribes two different possibilities. For small i?°^'°'s, the 
system is in a "smooth" phase where fluctuations in n° 's 
are suppressed and charges become mobile. The system 
is in a metallic phase. For large i?°^'°'s n°'s at different 
sites fluctuate violently (rough phase) and charge fluc- 
tuations are suppressed. The system is in the insulator 
phase [5]. 

Minimizing the free energy we obtain after some 
lengthy algebra the mean-field equations 



i?' 



J, 



cff ,o 
b 

MS 



Rab 



(7) 



= J a 



-(iA,.err> 



J; 



MS 



j: 



eff 



Jl 



off 



where P^°'{m 
{S{m — \n} — n: 
ger differences |A 
respectively. 



= {6{m - |A,nrl))s" and P^{m) = 
^|))sn are the probabilities that the inte- 



and 



trial' 



aa 1 2 \ 



(3Nd 



E 
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The phase diagram of the system is determined by solv- 
ing the above equations numerically. Notice that the 
superconducting transition given by = is deter- 
mined by S'triai only and is independent of Sj^iai ^ong 
as Pn°'{0) and /'^(O) are nonzero. Similarly, the metal to 
insulator transition is determined by S^^^.^^^ only (rough or 
smooth phase) when jf^ = 0. In Fig. 2 we present the re- 
sulting phase diagram for the symmetric case Jn = J22, 
Cii = C22 and i?ii = i?22 for two different values of J/. 
We note that a metallic phase is found in a narrow region 
of parameter space when J/ is small enough, contrary to 
the single-band case where no metallic phase is found. 
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Figure 2: Phase diagram of the two band system for different 
values of 011+022 versus _Ri2 for two values of Ji with jf^ 7^ 
in the upper panel and jf^ = in the lower panel. A metallic 
phase is found in a narrow region of parameter space when 
Jf = and 7?i2 < O.e/i/Qe^ 



where 1(2) = 2(1) and 



I 
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where 



2a/ 



abb = { Jbb + li^bb + ai2)\iOn\ ] l{k) + J/ 



1,2) and 



ai 



ai2 
2 



(8b) 



(9a) 



(9b) 



To understand the phase diagram, we observe first that 
the phase diagram is divided into two regimes, (i) J/^ 7^ 
(upper panel) and (ii)J/'^ = (lower panel). The effec- 
tive Josephson coupling between the two bands is nonzero 
in the first regime and it is easy to show from Eqs. 
([7]|9| that the two-band superconductor becomes effec- 
tively like a single band superconductor at low energy 
Correspondingly, Sll^^i becomes an 
Hi at temper- 
ature T ^ because of the pjf^^ ^ - Sin] — n|) term, 
i.e. 



w << Jf 



•^trial ^ o 2^ 



Qtt 



•off 



where aab = h/{ATre'^Rab)- The resistance RabS are given 
by^^QQ^ = 5(3i?aa^--RaB^)7 and R^^ = ^i2^^ + i(-f^i'r^ + 
R22^) where a°f, = h / (ATre'^ R^^) . This rather compli- 
cated form of resistance is a result of appearance of D^9 
terms in S'diss- j{k) = 4(sin^(fcx/2) -I- sin^(fcj,/2)) is the 
geometric factor of 2D square lattice. 



when superconductivity is destroyed, where Ct = C\\ -I- 
C22 and R^^ = R~{1 i?22^ '^Rxl ■ The system is 
at a roughening (insulator) phase when J^eff > Rci = 
0.6(/i/Qe^)[S] and a metallic phase exist only if — 
at a finite region of resistances Rd > ReS > RcS- For 
single-band superconductors, RcS = {h/Qe^) > i?c/[7], 
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and an intermediate metallic phase cannot exist in this 
case. 

We next consider regime (ii) where Jf^ — 0. First 

it is straightforward to show that Jj^ ^ as long as 

•^aa 0, indicating that the system behaves always like 

an effective one-band system at low enough energy in 

the superconducting state. The situation is different if 

superconductivity is destroyed. Substituting — 

into Eq. (8b), we obtain a self-consistent equation for 
Tofi — 



= k/|^„'(0)exp 



1 

'2(3Nd 



1 

'-.Jf + ai\uJnh{k) 



where aj 



ai 1^22 +012 (ail +022) 

I>ll+Ct22+4Q12 



(10) 



Equation ( 10 ) is solved numerically where we find 



that the equation has a non-zero solution only when 
|J/| > Jj(a7, Jii), which is a number depending on aj 
and roughly proportional to max( Jn, J22), the transi- 
tion from the Jj^ 7^ to Jj^ — state is a first order 



phase transition. The phase diagram determined by ( 10 ) 
is provided in the supplementary material. 

This interesting result suggests that although the su- 
perconducting state behaves always like an effective one- 
band superconductor at low enough energy, there ex- 
ists two kinds of non-superconducting states. The non- 
superconducting state is effectively one-band like when 
Jf ^ and two-band like when jf ^ 0. We find that 
an intermediate metallic phase may exist in the two-band 
like non-superconducting state. 

To see how this can occur we consider Eq. (8a I with 
Jf = 0. In this case we obtain 



(11) 



where 



h ( 1 



1 



corresponding to a single-band superconductor with ef- 
fective resistance = + {Raa + 2i?i2)~^, which 
is the effective resistance obtained from the resistance 
network model shown in Fig.(l). The SI transition is 
determined by Eq. (11) and S"^.^^^. To see the plausi- 



ble existence of metallic phase, we examine the limit 
R°2 ~ i?i2 — )■ 0. In this limit, a long-range order 
of H^^'s are built up in the winding number action 
•^triai because nj = nf^^Vi, and the system is always 
in the smooth phase. A metallic phase exists as long as 
i?eff ^ RiiR22/{Rii + R22) > h/Qe^ where ^ 0. 

The window for the existence of metallic phase nar- 
rowed down when R12 increases as shown in Fig. (2) 
lower panel. Notice that the winding number field 



is basically controlled by R12 when Rn and R22 are 
large, so for a metallic phase to occur, we generally re- 
quire R12 to be smaller than Q.&h/Qe^. For small i?i2, 
the superconductor-insulator transition is governed by 
ail + OL22- For the superconducting stiffness to vanish, 
we require an -I- a22 < 0.5. 

The metallic phase, if exists, is a new state of matter 
with incoherent bosons participating in low temperature 
transports. The state is described by a Ginsburg-Landau 
(GL) theory with vanishing phase-stiffnessfT^ . A prelim- 
inary analysis of the GL theory indicates that the system 
is a diamagnetic metal with unusual low-temperature 
magneto-transport behaviors jl2j. 

To conclude, we re-examine the problem of SI tran- 
sition in this paper for two-band superconductors, and 
raise again the question of plausible existence of metal- 
lic phase. Within a resistance-shunted Josephson net- 
work array model, we show that intermediate metallic 
phase between superconductor-insulator transition may 
exist for two-band superconducting films if the inter-band 
Josephson coupling J/ and inter-band dissipative resis- 
tance term i?i2 are small enough. Physically, the more 
complicated circuit network structure for two-band su- 
perconductors (Fig.l) gives rise to the possibility that 
the effective dissipation responsible for screening and 
quantum dissipation are coming from different resistance 
channels which is not possible for single-band supercon- 
ductors. With the recent advancements of research in 
Iron pnictide and other multi-band superconductors, we 
believe that this new metallic phase of matter may be 
reachable in the near future [TT]. 

This work is supported by HKRGC through grant 
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SUPPLEMENTARY MATERIALS 

Here we show how the dissipative term ([2| can be 
derived from a straightforward generaUzation of the dis- 
sipative term for one-band system to two-band system. 
We assume phenomenologically that a metaUic compo- 
nent exists in the system and the dissipative term can be 
derived from a Hamiltonian with tunneUing and capaci- 
tance energy between grains, 



(12) 



m,a a,b 



where m = L, R and a,b — 1,2 are the grain and band 
indices, respectively. 

Hn^ = E / rf^m*t"(^m)[C(-*V)]^^(x„) (13a) 
cr '' 

describes non-interacting electrons in grain m, band a 
where a is the spin index, and 

Hf^Y.J dxLdxRT''\xL,XR)¥j'{xL)i'U^R) + h.c. 

(13b) 

describes tunneling of electrons between grain L, band a 
and grain R, band b and 



1 



8C 



(13c) 



is the charging energy associated with charge imbalance 
between the grains where 

Ql^^Y. [ dx^¥f{x^)i'-M (14) 

is the total electric charge in grain to, band a. The cor- 
responding action at imaginary time is 

S [" dT{[ dxLKixL)dT-^l{xL) 



dxR^UxR)dT'^lixR)} + H 



(15) 



To derive Sdiss we first apply a Stratonovich-Hubbard 
transformation on Hn to obtain 



dx^^''„{x^){dT + e^^i-iV) 



1 + f-1)'''" 1 - f-l?"' 
+ ( \ ' )Va-a + ( Kr^)V,a]}KM 



- E] J '^^ y dx^dx} 

^b/_ \^l,a /_\,T,b , 



x[T'^\xL,XR,mXL)nixR) + C.C] 



-E/ ^-[^] 



(16) 



where m — L, R, = and sr = 1 and 1(2) = 2(1). 



Writing Vab = O^r- e\, where L{R) = R{L), the elec- 
tric potential VabS can be absorbed by a gauge transfor- 
mation 

Ki^ra^r) = e-''^"+i(''^+''-)*^(x„,T). (17) 
where the tunnelling term becomes 

St f '^'^ j dxLdxRiT'^^XL.XR) 

^gi[i(A„e°''+A„9'"')+i(A„ei2+A„e2i)] 

x^Uxl)'^''Axr) + c.c} (18) 

where A^6'''^ 6'£ ~ 6'^, {a,b = 1,2). To proceed fur- 
ther, we integrate out the fermionic fields and expand 
the tunnelling term to second order to obtain 



a,b 



' 1 

issJ 



(19) 



where 



oab 
diss 



Jei^^'i'/ 





dTl / dT2 

Jo 



{G2.(ri-T2)GL(r2"ri> 
+G''RAn-T,)GlAT2-n)e 



j(D„8'"'(Ti)-Z5„e'"'(r2)) 



} 



(20) 



where D^e}"^ = D^^f ^ = [A^Oj^ + A^Of] and D^Of" 



'^A„er + ^A^ef'']. 
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E 



sm^TrTr) 



(21) 

is the free electron Green's function at imaginary time. 
m — L,R and D'^^{Ep) is the density of states on the 
Fermi surface. Defining 



^sin(7rrT) ' 



h 



T 



-V 



4e^R°J Vsin(7rTr) / 
and put it back into pO] ) we obtain 

r/3 



(22) 



diss 



a,b 

xsin^[^-^"(-^);^^-^"(--)]} 



(23) 



which is the dissipation term we use in the main text. 

We attach here also the phase diagram determined by 
Eq. (10) with Jii = J22. The line separating the Jf^ — 
(7^)0 phases is a line of first order phase transition, we see 
aiia22+ai2(oii+a22) transition. 

Q:il+a22+4Qi2 



that Jill J I ^ a I 




Figure 3: Phase diagram for Jj^^ for different values of Jj 
versus a. A first order phase transition separates the Jj^^ = 
and Jj^^ =^ phases 



